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ABSTRACT 



In a five hour Ha exposure of the northwest region of the Coma cluster with the 2.1m telescope at San Pedro Martir (Mx), we 
discovered a 65 kpc cometary emission of ionized gas trailing behind the SBab galaxy NGC 4848. The tail points in the opposite 
direction of the cluster center, in the same direction where stripped HI had been detected in previous observations. The galaxy shows 
bright HII regions in an inner ring-like pattern, where the star formation takes place at the prodigious rate of ~ 8.9 M Q yr . From the 
morphologies of the galaxy and the trailing material, we infer that the galaxy is suffering from ram pressure due to its high velocity 
motion through the intergalactic medium. We estimate that ~ 4 x 10 9 M Q of gas is swept out from the galaxy forming the tail. Given 
the ambient conditions in the Coma cluster (p = 6.3 x 10~ 27 g cirr 3 ; a vel = 940 km s -1 ), simulations predict that the ram pressure 
mechanism is able to remove such an amount of gas in less than 200 Myr. This, combined with the geometry of the interaction, is 
indicative of radial infall into the cluster, leading to the conclusion that NGC 4848 has been caught during its first passage through 
the dense cluster environment. 

Key words. Galaxies: clusters: individual: Coma - Galaxies: individual: NGC 4848 - Galaxies: ISM - Galaxies: interactions 
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Yagi et al. (2010) reported a deep (4.5 hour integration) Ha 
survey covering the central 0.5 X 0.5 deg 2 of the Coma clus- 
ter (A1656) with the Suprime-Cam mounted on the Subaru tele- 
scope. Unexpectedly for an evolved cluster such as Coma, these 
authors found that almost every star-forming member has its own 
spectacular complex of diffuse, ionized, gaseous trails extend- 
ing dozens of kpc behind the optical extent of the galaxies, and 
sometimes harboring star-forming compact knots. They revealed 
14 such systems, including those previously reported by Yagi et 
al. (2007) and Yoshida et al. (2008) in the Coma cluster. Similar 
examples can also be found in A 1367 (Gavazzi et al. 2001), in 
Virgo (Yoshida et al. 2002, 2004; Kenney et al. 2008), and in 
A3627 (Sun et al. 2007). These features suggest that the galax- 
ies were recently captured by the cluster gravitational potential 
and are now infalling toward the cluster center (Yagi et al. 2010). 
Unfortunately, the field of view of the Suprime-Cam missed by 
less than 5 arcmin the position of NGC 4848 which is another 
obvious candidate for possible extended Ha emission. Bothum 
& Dressier (1986) had indeed listed NGC 4848 as one of the 
dozen unusually active galaxies found in the Coma cluster. 

NGC 4848 (CGCG 160-055; Zwicky et al. 1961-68) is a 
bright (M B =-20.5) SBab:edge-on (RC3, de Vaucouleurs et al. 
1991) galaxy that lies at the northwest (N-W) periphery of the 
X-ray emitting region in the Coma cluster. It has a vigorous star- 



* Observations taken at the observatory of San Pedro Martir 
(Baja California, Mexico), belonging to the Mexican Observatorio 
Astronomico Nacional. 



formation rate of ~ 9 M yr~' as derived from the Ha, ultravio- 
let (UV), far-infrared (FIR), and radio-continuum emission (see 
§0. Observations in the 21 cm line of HI (Gavazzi 1989; Bravo 
Alfaro et al. 2001) revealed a moderately deficient HI content 
(Defui = 0.46), displaced in the N-W direction, as opposed to 
its H2 content, which appears normal and centrally concentrated 
(Vollmer et al. 2001). The asymmetry in the HI distribution sug- 
gests that the galaxy is experiencing ram pressure (Gunn & Gott 
1972) owing to its high velocity motion through the intergalac- 
tic medium (IGM). This discrepancy is expected since unless 
ram pressure stripping is severe, only the atomic phase of the 
gas distributed at the galaxy periphery is removed, while the H2, 
bound deep within the galaxy potential well, is mostly unaffected 
by ram pressure (Combes et al. 1988; Kenney & Young 1989; 
Boselli et al. 2002; Fumagalli & Gavazzi 2008). 

Numerical hydrodynamical simulations of galaxies subject 
to ram pressure stripping in rich clusters (e.g. Kapferer et al. 
2009; Tonnesen & Bryan 2009, 2010, 2012; Ruszkowski et al. 
2012) reveal that in much less than 1 Gyr these galaxies lose all 
of their gas when the density of the IGM and the transit velocity 
are as high as in the Coma cluster. Consistent results are found 
both with and without magnetic fields. Extended gaseous tails 
form and the gas is shocked and heated by turbulence (Yoshida 
et al. 2004, 2012; Kenney et al. 2008), producing compact knots 
where radiative cooling takes place favoring the star formation. 

In 2000, we serendipitously discovered a low surface- 
brightness Ha emission trailing behind NGC 4848 (not reported 
by Iglesias-Paramo et al. 2002) in a one-hour exposure of the 
central lxl deg 2 of the Coma cluster with the Wide Field 
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Table 1. Log book of the imaging observations. 



Telescope 


Date 


CCD 


Pix 

(arcsec) 


Filter 

(A) 


Tint 

(sec) 


Nexp 


seeing 
(arcsec) 


INT 


20 Mar 1999 


4 x 2048 x 4100 EEV 


0.33 


B 


300 


2 


1.0 


INT 


28 Apr 2000 


4 x 2048 x 4100 EEV 


0.33 


6725 (80) 


1200 


3 


1.3 


INT 


28 Apr 2000 


4 x 2048 x 4100 EEV 


0.33 


(Gunn) f 


300 


3 


1.3 


TNG 


09 Feb 2001 


1024 x 1024 NICS 


0.25 


H 


60 


9 


0.8 


SPM 


Apr 2012 


1024 x 1024 EEV 


0.35 


6723 (80) 


600 


30 


1.4 


SPM 


Apr 2012 


1024 x 1024 EEV 


0.35 


(Gunn) r 


60 


26 


1.4 



Camera (WFC) at the Isaac Newton Telescope (INT, La Palma). 
However, this extended emission was only marginally detected 
and follow-up observations were required. Similar extended fea- 
tures were detected in deep GALEX images by Smith et al. 
(2010) showing several knots of recent star-formation along the 
tail. In 2012, we acquired additional five-hour observations with 
the San Pedro Martir (SPM) telescope using narrow-band Ha fil- 
ters. The resulting stacked six-hour exposure, which we present 
in this work, is sufficiently deep to allow a robust determina- 
tion of the flux in the tail. Throughout this paper, we assume 
H = 73 km s _1 Mpc~\ thus NGC 4848 is at the distance of 
95.5 Mpc, that of the Coma cluster. 



2. Observations 

We observed NGC 4848 in the B band using the WFC at the 
prime focus of the 2.5m INT in 1999 and in the H band using 
the Near Infrared Camera Spectrometer (NICS) mounted at the 
3.6m Telescopio Nazionale Galileo (TNG) in 2001. Both obser- 
vations (available through GoldMine, Gavazzi et al. 2003) were 
obtained under < 1.0 arcsec seeing (see Table [1]). We also ob- 
served the field centered on NGC 4848 in Ha using two tele- 
scopes: the INT, and the 2.1 m telescope at SPM. The observa- 
tions were performed through narrow-band filters centered at ~ 
6725 A, covering the redshifted Ha and [Nil] lines (ON-band 
image). The underlying continuum was measured through the 
broad-band r-Gunn filter (OFF-band image; see Table Q] for de- 
tails). The images were obtained in photometric conditions with 
seeing ranging from 1.3 to 1.4 arcsec. To minimize the impact 
of cosmic rays, we acquired a set of shorter exposures that were 
subsequently combined using a median filter. Images of the spec- 
trophotometric star Feige 34 were also collected for the photo- 
metric calibration. 

The individual images were bias-subtracted and flat-fielded 
using combinations of exposures of several empty fields at twi- 
light. A second-order flat-field correction was achieved using su- 
perflat images obtained by median-combining several science 
exposures that were selected because these are the ones contain- 
ing only point-like objects among all exposures taken during the 
2012 run. The images were aligned using field stars as refer- 
ence (using the IRAF task imalign) and fitted individually with 
third-order polynomials to correct for residual flat-field inhomo- 
geneities (using the IRAF task imsurfif). After background sub- 
traction and rebinning the two sets of images to a common pixel 
scale of 0.35 arcsec pixeL 1 (using the IRAF task imlintrari), we 
obtained the final combined image (using the IRAF task imcom- 
bine). The resulting stacked frames correspond in total to a 6 
hour and 41 min integration for the ON-band and OFF-band, re- 
spectively. 



Finally, we derive the continuum-subtracted Ha+[NII] im- 
age in the following way. First, the intensity in the combined 
OFF-band image was normalized to that of the combined ON- 
band one, using the flux ratio of several field stars. Then, since 
Spector et al. (2012) pointed out that the colors of foreground 
stars differing from those in the target galaxies can be an im- 
portant source of error when estimating the normalization coef- 
ficient, we applied a correction as follows. We evaluated that the 
average color of NGC 4848 is g - r ga i = 0.7 mag (see Table|2]i. 
This is only slightly different from the mean color of Milky Way 
stars at similar Galactic latitude g - r stars = 0.8 mag (Spector 
et al. 2012). From the linear fit shown in figure 2 of Spector et 
al. (2012), we infer that the normalization coefficient for NGC 
4848 is 0.98 times the one estimated from the foreground stars. 
In our case, this correction is only marginal. The combined NET- 
image was then obtained by subtracting the normalized OFF- 
band frame from the ON-band one. The resulting OFF-band and 
NET frames are shown in Fig.Q] The limiting Ha surface bright- 
ness was estimated using photometry of randomly sampled aper- 
tures of 8 arcsec in diameter. The resulting lcr limiting sensi- 
tivity is 1.0 x 10~ 17 erg cm" 2 sec _1 arcsec" 2 (at full resolution), 
which is a factor of four shallower than the superb Subaru data 
from Yagietal. (2010). 

3. The galaxy 

In spite of its SBab classification, the morphology of NGC 
4848 is rather irregular, being composed of a central nu- 
clear source (harboring an active galactic nucleus (AGN) 
[NII]k6584A/Ha=0.64, see Gavazzi et al. 2011) with red (g - 
i=\.5 mag) color. At a 5 - 10 arcsec projected radial distance 
from the nucleus, a bright, blue, and clumpy ring of active HII 
regions dominates the galaxy morphology. Further out, the light 
profile falls-off exponentially, as is typical of a disk. Given the 
irregular shape of the galaxy, a simple (bulge+disk) model can- 
not fit the measured light profile in any band. Therefore, a non- 
parametric classification method has to be used, such as the CAS 
system of Conselice (2003). 

Based on the Ha and r-band images obtained in 2012 and on 
B and H-band images, we derived the concentration, asymmetry 
and dumpiness (CAS) parameters of NGC 4848 (see Table 
We remind the reader that a theoretical r 1 ^ 4 law has C = 5.2, 
while an exponential disk has C - 2.7. Moreover, a symmetric 
and smooth galaxy has A — 0, S =0. Typical values for Sa-Sb 
galaxies in the r band are C = 3.9 + 0.5, A = 0.07 ± 0.04, and 
S = 0.08 ± 0.08 (Conselice 2003), which are consistent with the 
values of C = 3.0 + 0.5, A = 0.11 ±0.05, and S = 0.10 + 0.10 
found by Fossati et al. (in preparation) for a larger sample of 
360 galaxies in the Coma Supercluster in the r band. The same 
set of parameters were derived in the NET Ha frames, yielding 
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Fig. 1. Left panel: Grayscale representation of the Gaussian-smoothed (three-pixel-kernel) NET image of NGC 4848, with contours 
of the same image (contours are from 1 to 320 xlCT 17 erg cm -2 sec _1 arcsec~ 2 in 10 logarithmic intervals). The white cross marks 
the center of the molecular hydrogen given by Vollmer et al. (2001), while the white X marks the center of the HI given by Bravo 
Alfaro et al. (2001). Three knots of NET emission with associated stellar continuum are marked Kl-3. Right panel: Grayscale 
representation of the OFF image of NGC 4848, with contours of the Gaussian-smoothed (three-pixel-kernel) NET image (same as 
left panel). Note the angle between the major axis of the galaxy and the direction marked by the tail, which points toward the center 
of the Coma cluster. 



Table 2. Physical parameters of NGC 4848 that complement 
Table 1 of Vollmer et al. (2001). 



a (J2000) 


12"58'"05 s .6 


8 (J2000) 


28°14'34" 


Distance a 


95.5 Mpc 


U» 


15.82 mag 


«* 


14.54 mag 


r" 


13.81 mag 


f 


13.42 mag 


z" 


13.17 mag 


M Dya c 


3.8 x 10" M c 


M Sta / 


4.8 x 10 10 M c 


M m e 


2.2 x 10 9 M Q 


Mh/ 


4.1 x 10 9 M 


M ttail * 


3.6 x 10 9 M Q 


M B J 


5.8 x 10 10 M c 


M Dyn /M Bar 


6.5 


SFR* 


8.9 M yr-' 



Notes. ( '" Assuming a distance modulus in - M = 34.9 mag (Gavazzi 
et al. 2010). {b} Petrosian AB magnitudes from the SDSS DR7 database 
(Abazajian et al. 2009). Corrected for extinction in the Galaxy (Schlegel 
et al. 1998). <c) Assuming centrifugal equilibrium, V rot = 270 km s -1 
(Amram et al. 1992) and the B-band diameter «25 = 1.6 arcmin. (RC3, 
de Vaucouleurs et al. 1991) {d) Derived from the i magnitudes and g — i 
color using a modified Bell et al. (2003) recipe consistent with the mass 
determination of MPA-JHU. (e) Gavazzi (1989) (/) Vollmer et al. (2001) 
(s) This work (/,) M Bar = M StaI + M m + Mh, + M tai i 

C = 2.5 ± 0.9, A = 0.28 ± 0.11, and S = 0.65 + 0.26 for Sa- 
Sb. Therefore, NGC 4848 is significantly more asymmetric and 
clumpy and less concentrated than typical Sab galaxies, espe- 
cially in bands that more sensitively trace the presence of young 



stars, consistently suggesting that some morphological peculiar- 
ity has been recently (< 1 Gyr) induced by the interaction of the 
galaxy with its surrounding environment. 

Table 3. The concentration, asymmetry, and dumpiness (CAS) 
parameters of the galaxy NGC 4848 in four bands. 



Band C 


A 


S 


Ha 1.76 


1.05 


1.75 


B 1.83 


0.45 


1.06 


r 2.50 


0.30 


0.57 


H 3.06 


0.18 


0.35 



3.1. The SFR from the Ha luminosity 

To make a quantitative estimate of the Ha + [Nil] flux associ- 
ated with NGC 4848, we preferred to use the new SPM image 
alone to avoid uncertainties caused by to resampling that occurs 
when the image is stacked with the existing INT observation. 
The Ha + [Nil] flux within the optical extent of NGC 4848 was 
found to be logF(Ha + [Nil]) = -12.47 ± 0.03 erg cm~ 2 sec" 1 . 
Before converting this flux into a star-formation rate (SFR), the 
Ha flux must be corrected for Galactic extinction, contamination 
due to [NII]?Ji6548, 6584A, internal extinction, and dust extinc- 
tion, (see Boselli et al. 2009, Gavazzi et al. 2012). Given the high 
Galactic latitude of the Coma cluster, extinction due to the Milky 
Way (Schlegel et al. 1998) is negligible (the dust-corrected value 
islog J F(Ha+[NII]) MW = -12.46+0.03 ergcirT 2 sec" 1 ). The de- 
blending from [Nil] was computed using the drift-scan spectrum 



3 



M. Fossati et al.: Ha tail behind NGC 4848 




Fig. 2. Left panel: Grayscale representation of the inner structure of the NET image of NGC 4848 at full resolution, with contours of 
the OFF-band image. The white cross marks the center of the molecular hydrogen given by Vollmer et al. (2001), while the white X 
marks the center of the HI given by Bravo Alfaro et al. (2001). Right panel: High-contrast representation of the Gaussian-smoothed 
(three-pixel-kernel) NET image to highlight the emission in the tail, with contours of the OFF-band image. The black box marks 
the region where the flux in the tail was computed. 



of the galaxy taken from Gavazzi et al. (2004). The corrected 
flux became log F(Ha)faw-DB = -12.68 + 0.03 erg cm -2 sec -1 . 
The flux was further corrected for internal extinction using the 
Balmer decrement (as in Boselli et al. 2012, submitted), which 
was determined again with the drift-scan spectrum to yield 
log F(Hff)MW;DB;AA = -12.11 ± 0.03 erg cm -2 sec -1 . The fi- 
nal correction for dust absorption in HII regions assumes that 
only a fraction / of the Lyman continuum photons contribute 
to the ionization of the atomic hydrogen. Following Boselli et 
al. (2009) and assuming a Salpeter IMF, / = 0.77 we ob- 
tained log F(Ha)MW;DB;AA;EX = -11.99 ± 0.03 erg cirT 2 sec -1 . 
Conversely, the photon escape fraction was considered to be 
negligible (Boselli et al. 2009). The corrected Ha luminosity 
was thus Lfta = 1.1 + 0.2 x 10 42 erg sec -1 . According to 
Kennicutt (1998), the unobscured star-formation rate is given 
by SFR Ha [M yr 1 ] = 7.9 x 10 -42 x L Ha [erg sec' 1 ], thus 
S FRu a — 8.8 ± 0.2 M Q yr -1 . The bulk of the star formation takes 
place in bright HII regions distributed along the central ring-like 
pattern, as discussed above (see Fig.|2]Left). 



3.2. The SFR from the UV luminosity 



Using Kennicutt (1998) SFR UV [M G yr 1 ] = 1.4 x 10 -28 x 
Luv[erg sec 1 Hz -1 ], where Luv is obtained following Boselli 
(2011) from the GALEX far-ultraviolet (FUV) magnitude of 
16.29 mag (Martin et al. 2005). The UV flux is particularly af- 
fected by dust extinction and the correction (A(FUV)=2.01 mag) 
is obtained following Cortese et al. (2008) assuming that there 
is a balance between the absorbed UV photons and the total FIR 
(3-1 100 fim) energy. From the corrected UV magnitude, we ob- 
tained S FRw = 1 0.7 M yr - 1 . 



3.3. The SFR from the FIR luminosity 

According to Helou et al. (1985), the FIR 8-1000 fim luminosity 
is obtained by combining the 12, 25, 60, and 100 /urn flux from 
IRAS (Neugebauer et al. 1984). Since NGC 4848 is undetected 
in the IRAS 12 /mi band, we used the WISE (Wright et al. 2010) 
12 /urn flux. We obtained Fir = 1.51 x 10 -13 Wm -2 . Using the 
Kennicutt (1998) calibration SFR 1R [M G yr 1 ] = 4.5 x 10 -44 x 
LiR[erg sec -1 ], we obtained S FR\ R = 7.4 M yr -1 . 

3.4. The SFR from the radio continuum luminosity 

The estimate of the SFR from the radio continuum (RC) is 
based on the 1415 MHz flux (21.8 mJy) from NVSS (Condon et 
al. 2002) corrected for the nuclear contribution (7.1 mJy) from 
FIRST (Becker et al. 1994). Following the Bell et al. (2003) cal- 
ibration, SFR RC [M Q yr -1 ] = 5.52 x 10 -22 x L RC [W Hz -1 ], we 
obtained SFR RC = 8.8 M yr -1 . 

Assuming that the star formation is constant during the last 
10 7 (Ha)- 10 8 (UV, FIR, RC) yr (Iglesias-Paramo et al. 2004), the 
average star-formation rate becomes S FR — 8.9 ± 1.1 M yr -1 , 
corresponding to a specific star-formation rate (SSFR) of 1.9 x 
10 -10 yr -1 . This is more than a factor of ten higher than the 
SSFR of galaxies of similar mass and Hubble type. In conclu- 
sion, the galaxy appears to be undergoing a starburst phase pos- 
sibly triggered by the interaction with the cluster environment 
(as claimed by Bekki & Couch 2003). 

4. Ha tail parameters 

The flux estimate of the low-brightness tail requires a careful as- 
sessment of the quality of the flat-fielding and background sub- 
traction. For extended sources, the dominant source of error is 
associated with the background variations on scales similar to 
those of the source. We measured the residual background in 
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several 20 x 20 arcsec 2 regions (comparable with the extension 
of the trail) within the field and determined that the median of 
the residual background on this scale corresponds to 10 % of the 
background rms. From this figure, we estimated the uncertainty 
in the the Ha flux in the tail (see Gavazzi et al. 2012). 

At the distance of the Coma cluster, the total extent of the 
tail of 2.25 x 0.66 arcmin 2 corresponds to 62.5 x 18.5 kpc. 
To estimate the total flux in the tail, we integrated the counts 
above the 2-cr level in the full-resolution SPM image alone in- 
side a polygonal region of 2.25 x 0.66 arcmin 2 , excluding the 
galaxy itself and the compact knots (see below) (the polygo- 
nal region is sketched in Figure [2]). We obtained log F Vd n = 
-14.40 (-14.66) + 0.12 erg cm" 2 sec -1 . Hereafter, the value 
in parenthesis is obtained by correcting for [NII])Ji6548, 6584A 
contamination when a [NII]/Ha ratio of 0.62 is assumed, as for 
the main galaxy. Q No correction for internal extinction was ap- 
plied to the tail. 

The line intensity of the low-brightness tail was obtained us- 
ing 

7 H<* = , F ' ail =0.13 (0.07) Rayleigh ( 1 ) 

(1 Rayleigh = 10 6 /47r photons cirT 2 s" 1 sr" 1 ). Following Spitzer 
(1978), the Ha line radiation from an optically thin nebula at 
T = 10 4 K is customarily expressed in terms of the emis- 
sion measure (EM). Adopting an "effective recombination co- 
efficient" of a 32 = 1.17 x 10 _13 cm 3 s _1 (Osterbrock 1989), we 
derived lRyH<? = 2.77 EMhq;. 

In the case of NGC 4848, EM amounts to 0.36 (0.19) cm" 6 pc. 

Assuming that the tail has cylindrical symmetry, this implies 
a mass of ionized hydrogen of 5.0 x 10 9 (3.6 x 10 9 ) M , thus a 
plasma density of 4.3 x 10~ 3 (3.2 x 10~ 3 ) cm 4 . If the filling fac- 
tor is < 1, the values obtained should be taken as an upper limit 
to the mass and a lower limit to the density. The plasma in the 
tail must have originally consisted of neutral hydrogen, there- 
fore a rough estimate of the total gas loss from the galaxy can be 
derived using the HI deficiency parameter defined by Giovanelli 
& Haynes (1985). The galaxy contains 2.2 x 10 9 M (Gavazzi 
1989), corresponding to an HI deficiency parameter of 0.46, im- 
plying that 4.1 X 10 9 M of gas has been left behind in the tail. 
This estimate is consistent with the mass previously found in the 
ionized gas. 

We assumed a transit velocity through the cluster of ~ 
1330 km s _1 , as derived from V2<t,, c; , (where we also assumed 
that the cluster can be modeled as an isothermal sphere, as 
in Cayatte et al. 1994), where cr lW = 940 km s -1 for A 1656 
(Gavazzi et al. 2010). From the length of the tail, we then de- 
rived that the ionized material survived some 10 76 yr. The re- 
combination time r r = 1/N e ai2, is 10 78 (10 79 ) yr, i.e. consistent 
with the survival time. Although recent simulations favor a sce- 
nario where the gas is heated in situ to Ha emitting temperatures 
and the eddies along the tail contain sufficient turbulent energy 
to currently sustain its ionization (Tonnesen & Bryan 2010), it 
cannot be excluded that part of the gas in the tail was ionized by 
the galaxy HII regions or AGN. 

In addition to the ionized diffuse gas, the tail of NGC 4848 
harbors some compact star-forming regions (possibly HII re- 
gions) that we label Kl-3 in FigureQ] They are visible on the r 
band image and they coincide with compact features in the FUV 

1 The presence of a central AGN might lead to an overestimate of the 
[NII]/Hg' ratio. However, Yoshida et al. (2012) found that this ratio in 
the filaments associated with IC4040, another giant spiral galaxy in the 
Coma cluster, is that of a LINER-like spectrum. 



image by Smith et al. (2010). Their uncorrected Ha + [Nil] flux 
is 3.6x 10~ 16 , 2.1 x 10~ 16 , and 7.4 x 10~ 16 erg cm" 2 sec" 1 respec- 
tively. This corresponds to a luminosity of 4 X 10 38 , 2.3 x 10 38 , 
and 8 x 10 38 erg sec" 1 , which are typical of HII regions, and 
represent ~ 25% of the total luminosity in the tail. Similar "fire- 
balls" are found in the wakes of other stripped galaxies, such 
as GMP3016 and GMP4060 in the Coma cluster (Yagi et al. 
2010), or either VCC 1217 (Fumagalli et al. 2011) and VCC 
1249 (Arrigoni Battaia et al. 2012) in the Virgo cluster. These 
authors claim that these HII regions were recently born in situ 
out of stripped material, as predicted by some hydrodynamical 
simulations (e.g. Kapferer et al. 2009, Tonnesen & Bryan 2012). 

Boissier et al. (in preparation) claim the absence of in situ 
star formation behind galaxies in the Virgo cluster displaying 
HI (Chung et al. 2007) but no Ha tails. As explicitly mentioned 
by Tonnesen & Bryan (2010, 2012), the pressure by the IGM 
(e.g. density, velocity dispersion) is crucial for determining the 
formation of relevant eddies that might increase the IGM density 
enough to provide the cooling conditions driving the formation 
of molecular clouds, thus of stars. Such conditions are met in 
the Coma cluster and near the center of Virgo, but not further 
out where most of the Boissier et al. objects are located. 



5. Discussion 

There is unanimous consensus (Yagi et al. 2010, Yoshida et al. 
2012) that the morphological disturbances suffered by several 
late-type galaxies in the Coma cluster are caused by the dynami- 
cal interaction with the IGM, namely by the ram pressure mech- 
anism (Gunn & Gott 1972), including NGC 4848 (Gavazzi 1989, 
Vollmer et al. 2001). The question is how long this process has 
been acting for and whether the galaxy is infalling into the clus- 
ter for the first time or, as maintained by Vollmer et al. (2001), 
it entered the cluster environment about 1 Gyr ago and already 
crossed the cluster center 400 Myr ago. 

From the length of the tail, we inferred that the galaxy mo- 
tion takes place primarily in the plane of the sky, probably at 
V ~ 1330 km s _1 (y/2cr ve i). Unfortunately, we cannot directly 
measure this velocity since from the redshift we can infer only 
the component along the line of sight (V//). Subtracting from the 
heliocentric velocity of NGC 4848 (Zabludoff et al. 1993), the 
mean velocity of Coma cluster galaxies (Gavazzi et al. 2010), 
we obtained V// = 257 km s" 1 . 

There are two possibilities: that the orbit is circular, with its 
angular momentum in the plane of the sky, or that it is radial. 
In the first case, the ratio of the velocity along the line of sight 
to the tangential velocity V///V = 257/1330 = 0.19 gives the 
angle between V± and V. Assuming that the present projected 
radial distance of NGC 4848 from the cluster center is 0.75 Mpc, 
we conclude that the radius of the circular orbit is 4 Mpc, i.e. ex- 
ceeds the virial radius of Coma by almost a factor of two. We are 
much more likely to conclude, however, that the orbit is radial 
or at least highly eccentric, and that the present motion of NGC 
4848 is toward the center of the cluster at V ~ 1330 km s , as 
suggested by the direction of the Ha tail pointing in the opposite 
direction. 

More evidence of infall into the Coma cluster is provided 
by the velocity distribution of late-type galaxies (LTG) which 
in general hardly follow a Gaussian distribution but one that is 
skewed toward either lower or higher velocities than the overall 
mean velocity (Biviano et al. 1997, Boselli & Gavazzi 2006). 
The velocity dispersion profile of the LTGs in clusters is found 
to be consistent with orbits more radial than those of early-type 
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galaxies (ETG), providing a picture in which possibly all spirals 
have not yet crossed the virialized cluster core, and may even 
be on a first (infall) approach toward the central, high-density 
region. 

In the ram pressure scenario, the amount of stripped gas 
can be computed from the equilibrium between the gravita- 
tional binding energy and the dynamical pressure. The calcula- 
tion requires the density profile of the gas in the IGM. Following 
Cavaliere & Fusco-Femiano (1978), the density radial profile is 
well-approximated by an isothermal sphere 



PlGM - P() 



1 + 1^ 



-3/3/2 



(2) 
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where po represents the central density and r c the scale-length. 
Using the parameters estimates of Mohr et al. (1999) for the 
Coma cluster (all quantities having been recomputed for Ho = 
73 km s _1 Mpc~') namely po = 6.3 x 10~ 27 g cirT 3 , r c = 
0.26 Mpc, and /3 = 0.7, we computed that at the present posi- 
tion of NGC 4848 (0.75 Mpc from the center) the IGM density 
is po.75 = 6 x 10~ 28 g cirT 3 . 

The radius at which ram pressure becomes efficient can be 
estimated as (Domainko et al. 2006) 



^strip = 0.5R m ln 



( GM slR[ M m " 
v 2 PlCM 27TR 2 star R 2 m/ 



while the stripped mass is 



M strip = Mm 



^strip 
^HI 



+ 1 exp 



R 



strip 



R 



HI 



(3) 



(4) 



In this calculation, we assumed exponential profiles for both the 
stars and the interstellar gas, with /? star = 5.0 kpc (computed on 
the H band image) and Rm = 1.8x ^ star = 9.0 kpc for the HI 
disk (Boselli & Gavazzi 2006). This yielded /? stl -i p = 8.5 kpc and 
M strip = 4.7 x 10 9 M ~ 0.75 M H iori g - This is in remarkably 
good agreement with the mass of stripped gas (3.6 x 10 9 M ) 
computed in $4] which in turns is consistent with the missing 
mass of atomic hydrogen (4.1 x 10 9 M ) computed from the HI 
deficiency parameter. 

Using the ram pressure simulation by Kapferer et al. (2009) 
(the case with V re i = 1000 km s~') and assuming the IGM den- 
sity computed above (po.75 = 6 x 10~ 28 g cm -3 ), we derived 
that 65 % of the original HI content is stripped in about 200 
Myr. In spite of the different conditions found in the Virgo clus- 
ter with respect to the Coma Cluster, a similarly short timescale 
(~ 100 Myr) is reported by Boselli et al. (2006) for an almost 
complete ablation of the atomic gas from NGC 4569. This time 
is significantly shorter than the crossing time in the Coma cluster 
(1.6 X 10 9 yr, Boselli & Gavazzi 2006), supporting the conclu- 
sion that NGC 4848 is on its first passage through the cluster 
core. This time would be sufficient to remove most of its gas, in 
contradiction to the 1 Gyr timescale proposed by Vollmer et al. 
(2001). 
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